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The  detrended  data  base  has  been  analyzed  to  determine  the  probability 
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cumulative  probabilities  that  a  given  Q  is  exceeded  for  a  selected  duration 
(periods  from  0.  25  to  24  hr)  have  been  established  for  the  fully  detrended 
data  base.  To  interpret  the  results  in  their  effect  on  the  Over-the-Ilorizon 
Experimental  Radar  System,  the  distribution  of  probabilities  has  been 
modulated  with  the  dynamics  of  the  auroral  oval  bv  use  of  the  Starkov 
equation  (1969).  This  permits  the  assessments  of  the  probability  that  a  fixed 
location  in  the  radar  coverage  area  is  under  the  auroral  oval.  { 

For  determining  the  durations  of  auroral  disturbances  at  various 
corrected  geomagnetic  latitudes  within  any  test  area  unde>-  surveillance,  a 
number  of  tables  and  figures  are  presented  to  explain  the  procedure. 

For  the  convenience  of  a  user,  a  series  of  graphs  are  included  in  the 
report.  These  graphs  present  the  probability  of  auroral  disturbances  of 
various  durations  at  different  latitudes  in  the  Experimental  Radar  System 
(ERS)  test  area. 
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Magnetic  Disturbance  Statistics  From  a  Single  Station 
Q  Index  Applied  to  an  Actual  OTH~B  Radar  Situation 


1.  INTRODUCTION 

In  a  preliminary  study,  the  magnetic  activity  index  Q  (observed  from  a  single 
ground  station,  Sodankyla,  geographic  coordinates  67.4°N,  26.  corrected 
geomagnetic  coordinates  63.4°N,  108.  9°E)  was  used  to  determine  the  duration  of 
various  levels  of  magnetically  disturbed  periods.  It  was  noticed  that  the  Q  index 
measured  from  a  single  station,  which  moves  in  and  out  of  the  auroral  oval  due  to 
the  rotation  of  the  earth,  shows  a  strong  diurnal  variation.  This  report  attempts 
to  eliminate  the  diurnal  trend  fron.  the  data  and  presents  the  corrected  statistics 
of  the  durations  of  the  disturbed  periods.  This  analysis  is  based  on  the  Q  index 
data  for  the  year  1974  that  showed  the  largest  magnetic  disturbances  during  the 
20th  solar  cycle. 

In  the  second  half  of  this  report,  the  statistics  of  the  duration  of  the  disturbed 
periods  and  the  analytical  presentation  of  the  auroral  oval  are  combined  to  deter¬ 
mine  the  impact  of  the  aurora  on  the  geographical  area  of  the  OTH-B  Fxperimental 
Radar  Svstem  (KRSi. 


(Received  for  publication  18  December  1970) 
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2.  ANALYSIS  AMI  INT KR PR KT ATION 


To  determine  a  nmcpss  lor  the  plimination  of  thp  diurnal  variation  of  thr  Q 
index  from  a  single  station  (Sodankvla),  it  is  necessary  to  undprstand  the  cause 
and  the  extent  of  '.he  diurnal  variation.  For  this  purpose  the  decile  curves  for  the 
1!)74  Q-index  data  are  presented  in  Figure  1.  The  abscissa  shows  the  time  in  FT 
and  the  ordinate  shows  the  O  index.  The  ordinate  scales  of  the  curves  are  staggered 
to  avoid  the  mixing  of  the  decile  curves.  A  look  at  the  figure  reveals  that  all  the 
curves  show  a  systematic  minimum  between  06-10  FT  and  a  maximum  between 
20-02  FT.  Combining  the  auroral  dynamics  with  the  geographical  location  of 
Sodankvla.  results  in  the  maximum  distance  between  the  oval  and  the  station  at 
09  UT  and  the  minimum  distance  at  23  FT.  These  timings  match  with  the  timings 
of  the  minimum  and  of  the  maximum  in  the  diurnal  variation  of  the  Q  index.  Addi¬ 
tional  proof  is  the  behavior  of  the  top  percentile  (100  percent),  that  includes  the 
highest  values  of  the  Q  index.  The  "highest  values  refer  to  severe  magnetic  dis¬ 
turbances.  One  would  expect  these  to  occur  randomly  and  thus  be  uniformly  dis¬ 
tributed  over  the  time  scale.  As  shown  in  Figure  1,  the  Q  index  never  exceeded 
the  value  of  5  in  the  time  interval  0730-0913  UT,  whereas  it  attained  the  levels  8 
or  9  in  the  time  intervals  1900-0500  FT. 
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Thus,  even  during  the  most  severe  global  disturbances,  the  diurnal  oval  motion 
modulates  the  observed  level  of  disturbance  showing  that  Q  is  trulv  an  auroral 
parameter. 

The  distribution  of  the  observed  Q  values  for  197  4  is  similar  for  all  phases  of 
tile  solar  cycle.  The  fact  is  seen  in  the  100  percentile  curves  for  1965,  I960,  197  1 
and  (  1974  for  completeness)  again  1974  (figures  2a  to  2d'.  However  the  maximum 
level  of  Q-  9  is  never  observed  during  1965,  the  year  of  minimum  activity  in  the 
previous  (20th)  solar  cycle.  Also  the  highest  values  of  Q  observed  during  the 
quietest  time  of  the  day  (Q-3  at  0H1S  1'T1  are  lower  than  for  any  other  year  of  the 
solar  cycle.  The  most  enhanced  distribution  of  Q  is  found  in  1974,  about  three 
years  after  197  1,  the  vear  o!  ni.ixti’  urn  solar  activity. 

The  described  distribution  •>!  the  10H  per.  entile  curve  of  Q  with  a  minimum 
around  local  magnetic  midnight  was  first  indt.-e  1  bv  Martels  and  Fukushinia/  who 
initially  suggested  the  method  of  determining  the  Q  index  from  D,  anti  H  observa¬ 
tions  of  magnetic  components  from  turorai  stations. 

Thus,  one  can  visualize  a  correction  actoi  ,  c  m  nuent  upon  the  dc-i.  ,v.e  between 
the  locations  of  the  oval  and  the  stall  >o,  tor  ;he  e!it  fn-it’or-  of  the  d'jrr.a:  variation 
of  the  Q  index.  In  determining  the  fader,  several  ite  :s  e.-i.e  to  be  c<  ust  lered.  The 
Q  index  has  a  range  from  0  to  9  (also  1  and  F,  Ten  and  ldeven  rcspectivelv  ,  and 
an  unequal  distribution  with  60  percent  of  the  values  in  the  range  of  >  2.  Thus  a 
simple  arithmetical  average  would  not  be  appropriate.  To  overcome  the  difii  ultv 
of  eliminating  the  daily  trend  in  the  distribution  of  the  Q  index,  a  ratio  is  deter¬ 
mined  from  each  percentile  curve  by  dividing  the  highest  value  of  the  percentile  b  ■ 

(he  corresponding  15-tnin  value  of  that  percentile.  Percentile  curves  with  Q-  0  arc 
not  used.  Thus  referring  to  Figure  1,  of  the  ten  curves,  the  top  seven  with  nonzero 
Q  are  used.  The  effect  of  not  incorporating  the  lowest  three  decile  curves  into  ti  e 
deli  ending  process  is  small,  since  the  10  percentile  lias  no  diurnal  vari  item  and  tor 
the  20  an  )  30  percentiles  the  diurnal  variation  is  onlv  one  Q  level.  In  omparison 
the  diurnal  Q  level  ranges  for  the  higher  percentiles  are  A  Q-  1,  2,  3,  3,  4,  4  and 
5  levels  for  the  40  lo  100  percentile  curves  respective! v.  From  these  factors  an 
average  factor  for  each  15 -min  interval  is  computed  for  that  vear.  The  study  shows 
that  the  factors  for  each  vear  do  not  show  anv  obvious  dependence  on  die  phase  ol 
the  solar  cycle.  Therefore  for  eliminating  the  solar  cycle  dependence  if  anv,  factors 
are  averaged  from  the  data  for  the  four  years  1965,  1969,  197  1,  and  197  4.  The 
final  fa<  airs  are  plotted  in  Figure  3  and  are  .summarized  in  1  able  1. 

1.  Martels,  I.,  and  Fukushinia,  N.  (1956)  Kin  Q-index  fur  die  erd  magnetisi  he 
uktivitat  in  viertestundlidien  intervallen,  Ahhandl.  Akad.  Wiss.  (iotingen. 

Math.  -Phvs.  Klarse.  Sonderhefl -2. 
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detrending  factor 


Table  1.  Empirical  Factors  for  Detrending  the  Q  Index  Data  From  Sodankyla 


Hour  Quarterly  Intervals  of  the  Hour  Hour  Quarterly  Intervals  of  the  Hour 

(UT>  1234  1234 


0 

1.  0 

1.  1 

1.  1 

1.  1 

12 

1.  7 

1.  7 

1.  6 

1.  6 

1 

1.  1 

1.  2 

1.  3 

1.  3 

13 

1.  6 

1.  6 

1.  5 

1.  5 

2 

1.  4 

1.  4 

1.  5 

1.  6 

14 

1.  5 

1.  5 

1.  5 

1.  5 

3 

1.  6 

1.7 

1.  8 

1.  9 

15 

1.  4 

1.  4 

1.  4 

1.  4 

4 

1.  9 

1.  9 

2.  0 

2.  0 

16 

1.  4 

1.  4 

1.  4 

1.  4 

5 

2.  1 

2.  2 

2.  2 

2.  2 

17 

1.  4 

1.  4 

1.  3 

1.  3 

6 

2.  3 

2.  3 

2.  4 

2.  5 

18 

1.  3 

1.  3 

1.  3 

1.  3 

7 

2.  5 

2.  5 

2.  5 

2.  5 

19 

1.  3 

1.  3 

1.  3 

1.3 

8 

2.  6 

2.  6 

2.  6 

2.  5 

20 

1.  2 

1.  2 

1.  2 

1.  2 

9 

2.  5 

2.  5 

2.  4 

2.3 

21 

1.  1 

1.  1 

1.  1 

1.  1 

10 

2.  2 

2.  1 

2.  1 

2.  0 

22 

1.  0 

1.  0 

1.  0 

1.  0 

11 

2.  0 

1.  9 

1.  8 

I.  7 

23 

1.  0 

1.  0 

1.0 

1.0 

Multiply  the  observed  Q  from  Sodankyla  by  the  respective  factor  from 
the  table  for  detrending  the  diurnal  variation 


Figure  3.  Factors  for  Detrending  lit"  Observed  Q  Index  Data 


The  detrending  factors  for  each  15 -min  interval  are  plotted  in  Figure  3.  The 
abscissa  shows  the  time  of  observation  in  UT  and  the  ordinate  presents  the  factor 
for  correcting  the  observed  Q  value  for  the  diurnal  and  solar  cycle  dependence. 
Figure  3  shows  that  during  the  afternoon  and  nighttime,  the  factors  range  between 
1.  0  and  1.  5  but  the  factors  exceed  2.  0  for  a  6-hr  time  interval  centered  on  08  UT. 

In  the  following,  we  will  investigate  how  well  the  detrending  factors  in  Figure  3 
eliminate  the  diurnal  bias  from  the  observed  data  base.  The  factors  are  presented 
in  Table  1  so  as  to  be  easily  available  for  use. 

In  using  these  factors  \t  should  be  noted  that;  fat  the  observed  values  of  the 
Q  index  are  integers  ranging  front  0-9  (scale  extends  up  to  11);  who::  multiplied, 
some  values  may  exceed  the  range;  and  (b)  the  factor  will  not  change  the  observa¬ 
tions  with  Q=  0. 

This  second  shortcoming  can  be  partially  overcome  by  replacing  the  observe  i 
value  Q=0  by  Q=0.  25  and  rounding  the  value  aft  'r  the  multiplication  with  the  re  - 
spective  factor.  This  brings  in  a  dilemma.  Only  in  the  range  of  05-11  FT  ’factor 
i  2.0)  will  this  correction  have  an  effect  on  original  Q>  0  'h.-t  the  t  nr t'-m- tc  Q  in 
this  time  inverval  will  invariably  he  1.  However  the  replacement  <0-0.25'  will  not 
affect  the  original  values  of  Q=0  in  the  time  interval  I'-t/V.'T.  The  overall  effec 
in  such  replacement  Q=0.  25,  will  be  that  the  vet  y  mildly  disturbed  conditio' s.  with 
corre-tetl  0-1  in  the  interval  05-il  I  T  will  be  s'ightlv  overemphasized. 

In  view  of  t lie  dilemma  that  retaining  0-0  leads  to  an  underestimation  and  the 
replacement  of  Q=0  by  Q=0.  25  results  in  slight  overestimation,  and  in  both  cases 
to  a  small  distortion  of  the  Q  distribution  for  quiet  conditions,  0=0.25  has  been 
chosen. 

'i  lie  effect  of  applying  the  factors  to  correct  the  diurnal  variation  in  1974  is 
shown  in  Figure  4.  In  the  figure,  the  abscissa  is  the  time  in  UT  and  the  ordinate 
presents  the  Q  index.  The  continuous  curves  show  the  uncorrected  values.  The 
lower  curve  is  for  the  average  value  and  the  upper  curve  is  the  average  plus  one 
standard  deviation.  The  corresponding  detrended  curves  are  siiown  bv  the  dashed 
lines.  The  diurnal  variation  of  the  uncorreoted  average  Q  index  ranges  from  !.  2 
to  3.2.  The  averages  determined  from  the  adjusted  values  range  from  3.0  to  3.2. 
For  the  top  curve  (average  •  ct  )  the  range  for  uncorrected  Q  is  from  2.  2  to  5.  6. 

The  adjusted  (average  curve  *  o  )  ranges  from  5.3  to  6.2.  Thus  the  application  of 
the  orrection  .‘actors  significantly  reduces  the  range  of  the  diurnal  variation. 

The  reduction  in  the  diurnal  variation  should  also  be  evident  in  Hie  seasonal 
variation  of  Hie  <.)  index.  The  seasonal  variation  for  0!)  LIT  (minimum)  and  23  UT 
(maximum)  in  the  diurnal  variation  of  uncorreoted  values  and  detrended  values  is 
shown  for  the  year  ld/4  in  Figure  5. 
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Q  INDEX 


Figure  4.  Diurnal  Variation  of  the  Q  Index  From  Sodankyla  (Year  1974). 
Continuous  curves  are  for  uncorrected  data.  Dashed  curves  are  for  the 
detrended  data 


Figure  5.  Seasonal  Variation  of  the  Q  Index  From  Sodankyla  (Year  197“). 
Continuous  curves  are  for  the  uncorrected  data.  Dashed  curves  are  for 
the  detrended  data 


In  Figure  5,  the  lower  continuous  curve  shows  the  seasonal  variation  of  the 
observed  uncorrected  Q  values  for  09  UT.  The  upper  continuous  curve  shows  the 
variation  for  23  UT.  The  difference  between  the  curves  is  minimum.  AQ*1.5  in 
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July  and  maximum,  AQ=2.  6  in  April.  When  the  observed  Q  values  are  adjusted  by 
using  the  factors  (from  Table  1),  the  curve  for  23  hr  UT  is  unaffected,  since  this 
hour  provided  the  reference  values  for  the  detrending  factors.  The  adjusted  curve 
of  0!)  in,  in  dashes,  lies  now  close  to  the  curve  of  23  UT  with  largest  AQ=1.3  in 
the  month  of  July.  For  all  other  months  AQ  is  less  than  0.  5.  Again  the  factors 
reduce  the  diurnal  variation  in  the  observed  uncorrected  data  for  all  the  seasons. 

Due  to  the  systematic  minimum  in  the  diurnal  variation  in  the  observed  values 
of  the  Q  index  (uncorrected  data),  the  duration  of  a  disturbance,  exceeding  a  given 
Q  value  based  on  the  original  data,  will  be  underestimated.  In  fact  this  is  seen  by 
dividing  the  data  in  two  12-hr  intervals,  one  around  the  period  of  the  minimum 
(0130-  1330  UT)  and  the  other  around  the  period  of  the  maximum  (1330-0130  UT)  o' 
the  diurnal  variation.  These  results  are  presented  in  Table  2,  from  the  data  of  107  1. 

Tables  2a  to  2f  present  the  probabilities  for  the  magnetic  activity  to  exceed  the 
levels  of  Q  from  0-8.  In  each  table  the  top  line  describes  the  selec  tion  of  the  flat.) 
base,  nighttime  only,  daytime  only,  or  all  the  data,  fables  2a  to  2c  are  derived 
front  the  uncorrected  data  and  Tables  Je  to  2f  are  from  the  data  detrended  by  the 
factors  in  Table  1.  The  column  on  the  left-hand  side  presents  the  level  of  the 
disturbance  equal  to  or  greater  than  the  given  Q  value.  There  are  25  columns,  each 
refers  to  the  minimum  duration  for  the  disturbance  for  the  Q  level  defined  in  the 
left-most  column.  Thus  the  first  column  presents  all  durations  for  a  given  Q  ievet, 
whereas  the  last,  column  on  the  right-hand  side  is  for  durations  exceeding  one  dev 
(2-1  hours),  'file  probability  for  the  occurrence  js  presented  as  a  percent  of  (be 
year.  The  top  line  for  each  section  is  100  percent  for  Q  -  0.  ’I  his  line,  though  an 
obvious  result,  is  added  for  the  sake  of  consistency,  and  is  needed  in  future  com¬ 
putations  . 

in  determining  durations  from  day  or  nighttime  data  the  succeeding  selected 
(dav/mght)  periods  are  treated  as  continuous.  This  continuity  enables  the  deter¬ 
mination  of  the  probability  of  durations  longer  than  ’  2  hr  (up  to  24  hr). 

in  general  in  a  given  table  the  probability  decreases  as  '.lie  duration  of  disturb¬ 
ance  is  increased  (along  a  given  row).  Also  for  a  given  duration  the  probability 
decreases  as  the  (.)  level  of  the  disturbance  is  increased  (down  a  given  column). 

A  comparison  of  probabilities  in  Tables  2a  and  2b  shows  that  the  corresponding 
probabilities  at  nighttime  are  larger  than  those  at  davtime.  The  nighttime  proh  - 
abilities  from  raw  data  (Table  2a),  show  no  significant  differences,  from  those  from 
the  tetreuded  data  (Table  2  1)  as  anticipated  in  the  study  of  'I  able  1.  The  probabilities 
for  the  davtime  and  those  for  the  combined  dav  and  night  (all  data  base)  show  larger 
changes  Cl  ables  2h  and  2c  and  2c  and  2D.  Another  interesting  feature  is  that  the 
con  e,  (c  l  (ail  data)  probabilities  (Table  2c)  are  about  double  of  those  uncorrected 
nighttime  probabilities  (Table  2a).  Thus,  on  the  whole,  not  correcting  for  'lie 
systematic  diurnal  variation  of  the  observations  results  in  a  severe  underestimation 
of  t lie  real  probability. 
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Another  approach  of  overcoming  this  difficulty  is  to  supplement  observations 
from  one  station  by  those  from  another  station  with  a  similar  latitude  but  with  a 
longitude  separation  of  approximately  180°.  The  Observatory  at  College,  Alaska 
(geographic  64.9°N,  212.  2°E)  satisfies  this  requirement  for  the  Observatory  at 
Sodankyla  (geographic  67.4°N,  26,  6°E). 

The  Q  values  from  College,  Alaska  are  not  routinely  available.  Therefore  an 
effort  is  in  progress  to  determine  the  Q  index  from  the  College  observations  of  the 
D  and  H  components.  Both  approaches,  detrending  the  single  station  data  and 
complimenting  single  station  data  from  a  second  station,  will  be  followed  for  over¬ 
coming  the  systematic  diurnal  variation  of  Q  from  a  single  station. 


3.  DETERMINATION  OF  THE  IMPACT  OF  MAGNETIC  DISTURBANCES  ON 

THE  REGION  COVERED  IN  THE  TESTS  OF  THE  OTH  RADAR  SYSTEM 

The  irregular  features  of  the  auroral  and  subauroral  ionosphere,  such  as  the 
auroral  F-layer,  auroral  E  the  highly  irregular  F-layer  (FL1Z  or  E-Laver 
Irregularity  Zone),  auroral  absorption,  the  mid-latitude  F-layer  trough,  and  the 
strong  gradients  at  the  interfaces  of  the  various  ionospheric  structures,  are  all 
organized  by  a  :ow  well  established  frame  of  reference;  The  auroral  oval.  ^  The 
dynamic  behavior  of  the  auroral  oval,  controlled  by  the  interplanetary  magnetic 
field  and  geomagnetic  activity,  is  reflected  in  associated  dynamic  changes  of  the 
location  and  intensity  of  the  various  ionospheric  features  mentioned  above. 

The  Over -the -Horizon  Backscatter  Experimental  Radar  System  (OTH-B  ERS) 
presently  under  construction  in  Maine  will  operate  in  an  area  covering  a  part  of 
Greenland,  bounded  in  the  east  by  Iceland,  and  extending  over  a  large  part  of  the 
North  Atlantic  Ocean. 

Fvery  night  the  auroral  oval  and  thus  the  auroral  ionosphere  will  affect  at  least 
half  of  the  F  RS  coverage  area.  The  coverage  area  is  shown  m  Figure  6.  super¬ 
imposed  on  a  section  of  the  northern  hemisphere,  which  is  plotted  in  corrected 
geomagnetic  coordinates.  This  coordinate  system  is  especially  useful  for  the 
presentation  of  auroral  phenomena  such  as  the  auroral  oval. 

Two  positions  of  the  statistical  auroral  oval  for  Q-  8  are  shown  in  Figure  6. 

At  02  FT  (23  LT/21  CGT  in  the  center  of  the  ERS  coverage  area)  the  oval  extends 
the  furthest  south  (64’  Corrected  Geomagnetic  Latitude)  and  as  shown  bv  the  location 
of  equatorward  border  (bold  continuous  curve)  it  covers  the  top  (poleward)  one -half 
area  of  the  radar  coverage.  At  14  FT  (11  I.T/00  CGT  in  the  center  of  the  coverage 
area)  the  night  sector  of  ttie  oval  is  in  the  opposite  quadrant  and  the  equatorward 
border  of  the  day  sector  shown  by  hold  dashed  curve  now  covers  onlv  a  small  section 

2^  f'el'lstem,  Y .  I.  ( 10  63 )  On  the  morphoh >gv  of  auroral  and  magnetic  disturbances 
at  high  latitudes,  Geomag.  Aeron.  3;  183. 

- - — _ _ 2 _ , _ wW 
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of  OTH-B  FRS  surveillance  area  (72.5°  CG1.I.  Referring  to  Figure  6,  it  is  seen 
that  for  any  given  level  of  magnetic  activity  (in  this  case  Q=3),  the  obscuration  of 
the  FRS-test  area  depends  on  time  due  to  the  diurnal  motion  of  the  oval. 


Figure  6.  Coverage  Area  of  the  OTII-B  FRS  System  in 
Geomagnetic  Coordinates.  Note  the  change  in  the  position 
of  the  (.}-  ;l  oval  at  Oft  (dashed  curve)  and  2  1  CG’i 
(continuous  curve) 


li» 


A  given  ground  station  at  auroral  Latitudes,  as  well  as  the  coverage  area  of 
the  ERS  as  a  whole,  moves  underneath  and  out  again  from  under  the  oval,  each 
24-hr  period.  Due  to  this  diurnal  motion,  the  auroral  oval  will  spread  over  the 
varying  sections  of  the  coverage  area  even  for  a  steady  Q  index  (as  seen  in  Q=3 
in  Figure  6).  Additional  variation  of  the  area  under  the  auroral  oval  influence, 
results  from  the  dependence  of  the  oval  diameter  on  Q.  These  two  variations  can 
be  determined  in  two  steps:  The  variation  due  to  (1)  the  daily  motion  of  the  oval  for 
each  Q  combined  with,  (2)  the  expected  variation  of  the  Q  index  based  on  measured 
probabilities. 

First  we  will  determine  the  duration  for  which  the  equatorward  boundary  of 
the  oval  will  stay  equatorward  of  a  given  CG  latitude  for  a  given  Q  level  of  magnetic 
activity.  Since  the  oval  exhibits  a  24-hr  periodicity  the  duration  is  expressed  as  a 


percent  of  a  24-hr  period. 
Starkov  which  describes  6 


For  this  purpose  we  will  use  the  analytical  equation  of 
,  the  location  of  the  equatorward  boundary  of  the  oval 


as  a  function  of  t,  the  hour  angle,  and  Q  the  level  of  magnetic  activity.  This  equa¬ 


tion  is  based  on  the  study  of  discrete  auroral  arcs.  The  equation  is  given  by 


e  =  18  +  0.  9  Q  +  5.  1  cos  (t-  12°) 
eq 


0  is  the  oval  boundary, 

Q  is  the  level  of  magnetic  activity  on  the  midnight  meridian  at  65°  CGL,  and 
t  is  the  hour  angle  measured  in  degrees  from  the  midnight  meridian  eastward 
towards  the  morning  sector.  The  equation  provides  an  instantaneous  position  of  the 
oval,  depending  on  the  value  of  the  Q  index.  The  percent  of  each  day,  for  which  a 
point  with  a  specified  CG  latitude  is  under  the  auroral  oval  has  been  determined 
from  Eq.  (1).  This  is  done  by  computing  the  two  values  of  angle  t  for  a  given  Q  and 
a  given  CG  latitude  from  Eq.  (1).  The  difference  between  the  angles  divided  by  360° 
provides  the  fractional  period  of  a  day  for  the  given  Q  and  A.  These  are  presented 
as  a  function  of  Q  in  Table  3. 

The  left-hand  side  column  presents  the  CG  latitude  of  the  equatorward  boundary 
of  the  auroral  oval.  The  Q  index  levels  are  shown  along  the  top  row.  The  fractions 
of  24-hr  periods  for  given  Q  and  A  are  presented  as  percentages  along  respective 
rows  and  columns.  For  example,  if  the  magnetic  activity  is  continuously  at  level 
Q=9,  the  equatorward  boundary  would  be  at  or  lower  than  63°  CGL  for  0.  5  for  each 
day,  that  is,  for  12  hr  of  each  24-hr  period.  The  latitudes  less  than  57°  are  not 
shown  since  the  statistical  oval  does  not  reach  such  low  latitudes.  In  actuality,  on 
rare  occasions  the  aurora  has  been  observed  at  latitudes  much  further  south  than  the 
57°  CGL  predicted  from  Table  3. 

Sta rkov,  G.  V.  (1969)  Analytical  representation  of  the  equatorial  boundary  of  the 
oval  auroral  zone.  Geomag.  Aeron.  (English  Edition)^L6 14. 
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Table  3.  Percent  of  a  Day  That  a  Point  With  a  Specified  CG  Latitude  A  is  Poleward 
of  the  Equatorial  Oval  Boundary,  as  a  Function  of  Q 


■ 

Q  Index 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

57 

6 

58 

8 

21 

59 

13 

21 

29 

60 

13 

23 

29 

38 

61 

15 

25 

31 

38 

44 

62 

17 

25 

33 

38 

46 

50 

63 

17 

25 

33 

40 

46 

50 

56 

64 

17 

27 

33 

40 

46 

50 

58 

63 

65 

21 

29 

33 

42 

46 

52 

58 

63 

7  1 

66 

6 

21 

29 

35 

42 

46 

52 

58 

65 

7  1 

79 

67 

21 

29 

35 

42 

48 

54 

58 

65 

7  1 

79 

94 

68 

29 

38 

42 

48 

54 

58 

67 

7  1 

79 

100 

100 

69 

38 

42 

50 

54 

60 

67 

73 

83 

100 

100 

100 

70 

44 

50 

54 

60 

67 

75 

83 

100 

100 

100 

100 

7  1 

50 

54 

63 

67 

75 

83 

100 

100 

100 

100 

100 

72 

56 

63 

69 

75 

85 

100 

100 

100 

100 

100 

100 

73 

63 

7  1 

77 

88 

100 

100 

100 

100 

100 

100 

100 

74 

7  1 

79 

88 

100 

100 

100 

100 

100 

100 

100 

100 

75 

79 

92 

100 

100 

100 

100 

100 

100 

100 

100 

100 

76 

94 

100 

100 

100 

100 

100 

100 

100 

100 

100 

100 

77 

100 

100 

100 

100 

100 

100 

100 

100 

100 

100 

Table  4a  lists  the  probability  that  the  observed  Q  would  lie  in  a  given  interval. 
Table  4b  lists  the  cumulative  probability  that  the  observed  Q  would  be  2  the  given 
value.  For  the  sake  of  completeness  the  values  are  listed  for  the  years  1965,  1969, 
1971,  and  1974  for  the  data  from  Sodankyla.  The  letter  (u)  refers  to  the  probabilities 
from  the  uncorrected  data  and  (d)  refers  to  those  from  the  detrended  data.  It  is 
seen  that  1974  has  the  highest  values  for  the  probability  (Table  4b). 

The  product  of  percent  period  from  Table  3  and  the  probability  in  Table  4b 
gives  the  average  percent  of  time  that  the  auroral  oval  is  equatorward  of  a  given 
latitude  in  the  coverage  area  (averaged  over  the  year  1974)  for  a  Q  value  greater 
than  the  selected  number. 

Based  on  Tables  3  and  4b  the  probable  average  annual  duration  of  coverage  by 
the  aurora  in  percent  of  time  in  the  OTH-B  EHS  coverage  area  is  shown  in  Figure  7. 
Starting  at  the  lefthand  top  corner  of  the  coverage  area,  a  very  small  region  shaded 
is  affected  all  the  time.  Approximately  one-quarter  of  the  area  is  affected  half  of 
the  time,  whereas  the  area  shaded  with  lines,  which  forms  one-quarter  of  the 
sector  in  the  right-hand  bottom  corner,  is  statistically  free  of  any  auroral  disturb¬ 
ances  at  all  times. 
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Table  4.  Probability  for  the  Distribution  of  the  Q  Index  for  Various 
Years,  (a)  Uncorrected  and  detrended,  (b)  Cumulative  Uncorrected 
and  detrended 


(a)  Probability  (%)  for  Q  to  be  in  a  Given  Interval 


Q 

Index 

(u) 

1965 

<d) 

(u) 

Year 

1969 

<d) 

(u) 

1971 

(d) 

(u) 

1974 

(d) 

0 

62 

43 

50 

36 

45 

30 

23 

15 

1 

19 

28 

19 

24 

22 

25 

22 

17 

2 

9 

11 

12 

12 

14 

13 

18 

13 

3 

5 

7 

8 

10 

9 

10 

15 

13 

4 

2 

4 

4 

7 

4 

7 

9 

12 

5 

1 

3 

3 

5 

3 

6 

6 

10 

6 

1 

1 

2 

3 

2 

3 

4 

6 

7 

1 

1- 

2 

1 

3 

2 

6 

8 

1 

2 

1 

4 

9 

2 

10 

1 

(b)  Cumulative  Probability  for  Q  to  be  i  the  Given 
Interval 


0 

100 

100 

100 

100 

100 

100 

100 

100 

1 

38 

57 

50 

65 

55 

70 

77 

85 

19 

28 

31 

41 

34 

45 

55 

68 

10 

18 

18 

30 

19 

32 

37 

55 

5 

11 

10 

20 

10 

22 

22 

42 

5 

3 

6 

6 

13 

6 

15 

13 

30 

6 

1 

3 

3 

8 

3 

9 

7 

20 

7 

2 

1 

5 

1 

6 

3 

14 

8 

1 

3 

3 

1 

8 

9 

1 

1 

4 

L0 

1 

1 

2 

u  -  Uncorrected.  d  -  Detrended 


In  Figures  6  and  7,  the  curves  have  similar  shapes.  In  Figure  6,  the  curves 
present  instantaneous  positions  of  the  oval  (as  an  example),  whereas  in  Figure  7 
the  curves  present  the  average  percentages  of  the  auroral  disturbances,  affecting 
different  regions  in  the  test  area.  As  these  two  figures  present  entirely  different 
aspects  of  the  effects  of  Q  on  the  test  area,  the  incidental  similarity  and/or  a 
match  in  the  shapes  of  the  curves  cannot  be  used  for  determining  the  average 
distribution  of  the  Q  index. 


F  igure  7.  Percent  Obscuration  of  the  OTH-B  EKS  Coverage 
Due  to  Auroral  Disturbances  in  1974  as  Determined  F  rom 
Sodankyla  Detrended  Q  Index 


In  Table  2  we  considered  the  durations  of  disturbances  which  depended  onlv  on 
the  level  of  the  O  index.  Table  3  gives  the  duration  (in  percent  of  a  davl  that  a 
given  point  is  under  the  auroral  oval,  as  a  function  of  Q.  In  order  to  determine 
the  probability  of  the  duration  of  a  disturbance  lasting  at  a  given  latitude  for  a  given 
length  of  time,  these  two  factors  have  to  be  combined.  This  global  probability 
independent  of  longitude  can  be  computed  from  the  equation 
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(2) 


PA,d  zP\.Q  <PQ..d'P  Q.+  1,  d) 


where 


^  is  the  probability  for  duration  d  at  geomagnetic  latitude  A, 

q  is  the  probability  for  the  oval  to  be  at  latitude  for  given  Q.  (from 
Table  3),  1 


Q.,c 

l 

Table  2f), 


is  the  probability  for  the  duration  d  for  the  given  Q  index  Q.  (from 


A  is  the  CG  latitude  of  a  given  location, 

Q  is  the  ith  detrended  Q  index, 

l  is  used  for  summation  over  all  the  values  of  Q  (for  Q  2  0,  F  j  =  1.0). 

These  probabilities  for  various  lengths  of  disturbances  occurring  at  different 
latitudes  are  listed  in  Table  5.  The  lefthand  column  lists  the  CG  latitude.  The 
top  row  lists  the  periods  for  the  duration  of  the  disturbances.  For  example  at  ~>0" 
CGL  durations  less  than  I  hr  have  occurrence  probability  of  1  percent  (0.  01)  where¬ 
as  at  77°  CGI,  the  disturbance  is  present  all  the  time  [1.0  (100  percent)!.  This 
computation  does  not  take  into  account  the  longitude  of  any  given  location.  Essen¬ 
tially  the  computation  converts  the  (Q,  d)  domain  of  Table  2  to  the  (A,  d)  domain  on 
the  basis  of  the  Starkov  equation.  In  computing  the  probabilities  with  the  use  of 
Eq.  (2),  the  respective  probabilities  from  Tables  2f  and  i  have  to  be  used  as  frac¬ 
tions  («  D  instead  of  percentages  presented  in  the  tables,  lor  example  to  compute 
the  probability  that  a  point  with  A  -  70°  CGI.  Is  under  a  disturbance  of  a  minimum 
duration  of  12  hr,  Eq.  (2)  uses  data  for  d  irations  ol  12  hr  (13th  column)  from 
Table  2f  and  data  for  A  -  70  (  14th  row)  from  I  .hie  2.  The  final  sum  is  presented 

as  percent  in  Table  a  along  A  =  70°  (12th  row)  and  duration  of  12  hr  (13th  column). 
‘The  result:  A  probability  of  a  1  percent  shows  that  the  equate rward  boundary  of  the 
oval  would  be  at  or  below  a  t'G  latitude  of  70°  with  a  continuous  duration  of  at  least 
12  hr  for  ’>  1  percent  of  the  vear.  This  probability  is  independent  of  the  longitude 
of  a  given  station.  For  determining  probabilities  for  a  test  coverage  area,  a  fixed 
geographic  longitude  of  the  location  on  the  earth,  has  to  he  taken  into  account. 
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In  deriving  probabilities  listed  in  Table  5  the  durations  at  various  latitudes 
from  Table  3  are  used.  It  is  recalled  that  these  durations  refer  to  the  maximum 
available  period  at  a  given  latitude  independent  of  longitude.  In  short,  it  refers 
to  the  farthest  extent  of  the  oval  without  restricting  or  limiting  the  location  to  a 
selected  longitude.  For  example,  referring  to  Eq.  (1)  for  the  Starkov  auroral 
circle,  the  oval  is  consistently  at  lower  latitudes  for  longitudes  of  the  night  sector 
and  is  consistently  at  higher  latitudes  for  the  longitudes  which  are  in  the  day  sector. 
Therefore,  the  probabilities  in  Table  5  really  refer  to  the  disturbances  affecting 
the  locations  on  a  global  scale.  In  practice,  one  is  restricted  to  a  limited  region  of 
a  selected  longitude  sector.  That  the  probabilities  in  such  a  case  are  much  smaller 
(especially  when  longer  durations  are  considered)  is  seen  in  the  following  approach. 

To  determine  the  probabilities  of  disturbance  durations  for  points  along  a 
selected  meridian,  a  direct  approach  was  used.  With  the  use  of  detrended  Q  index 
data  of  1974,  the  latitudinal  position  of  the  equatorward  oval  boundary  for  each  Q 
value  at  the  Narssarssuaq  meridian  (NQ  in  Figure  6)  was  computed,  using  the 
Starkov  equation.  From  this  table  of  boundary  locations,  the  probability  distribu¬ 
tions  of  durations  of  disturbances  as  a  function  of  CG  latitude  for  various  time  inter¬ 
vals  are  computed.  These  are  summarized  in  Table  6.  The  format  of  Table  6  is 
identical  with  that  of  Table  5. 

A  comparison  of  the  corresponding  columns  from  Tables  5  and  6  shows  that  the 
respective  magnitudes  of  probabilities  for  a  global  base  are  larger  than  those  for 
a  single  station  like  Narssarssuaq.  The  differences  are  smaller  for  shorter  dura¬ 
tions  in  the  columns  on  the  left-hand  side  of  the  tables  (<  5  percent  up  to  4  hours'. 
These  differences  increase  progressively  from  the  left-hand  side  toward  the  right- 
hand  side  of  the  columns.  The  differences  are  very  significant  for  larger  durations 
of  disturbances  (right-hand  side  columns).  For  example  at  75°  CGI.  the  probabilities 
for  the  disturbances  with  durations  of  3,  and  34  hr,  on  a  global  base  are  94  percent, 
and  37  percent.  The  respective  values  for  Narssarssuaq  longitude  are  91  percent, 
and  96  percent  respectively. 

1  or  determining  the  extent  to  which  the  probabilities  in  Table  G  depend  upon 
longitude,  similar  tables  were  computed  for  the  longitudes  of  stations  CAP  TOBIN 
and  C\KIBOU  (TZ  and  CM  in  Figure  G)  which  lie  on  the  end  longitudes  defining  the 
sector  of  the  EHS  test  area.  The  difference  in  longitude  between  these  two  stations 
is  equivalent  to  a  time  difference  of  5  hours.  The  comparison  of  corresponding 
values  wdth  those  in  Table  G  produced  an  excellent  agreement  (within  3  percent)  for 
95  percent  of  these  values.  Thus  one  can  conclude  that  the  values  in  Table  6  arc 
reasonably  independent  of  the  longitude  of  the  location  for  durations  up  to  21  hours. 


For  understanding  the  relation  between  the  corresponding  values  of  probabilities 
listed  in  Tables  5  and  6  one  has  to  consider  several  factors.  One  of  the  considera¬ 
tions  is  the  variation  of  the  (detrended)  Q  index  with  various  time-spans.  The  study 
of  the  1974  detrended  Q  index  data  from  Sodankvla  shows  that  on  the  average  the  Q 
index  changes  by  one  unit  in  the  interval  of  an  hour.  For  all  other  intervals  in  the 
range  of  2  to  24  hr,  the  Q  index  on  the  average  changes  by  two  units.  Referring  to 
the  Q  index  term  in  the  Starkov  Eq.  (1),  it  is  seen  that  a  unit  change  in  the  Q  index 
is  equal  approximately  to  1°  change  in  the  latitudinal  position  of  the  aurora.  Another 
factor  is  the  diurnal  motion  of  the  oval.  Again  from  the  Starkov  equation  one  finds 
that  at  quiet  magnetic  conditions,  that  is,  Q=0  the  Oval  latitude  changes  from  77° 
around  local  magnetic  midnight  to  67°  around  local  magnetic  noon.  Superposed  on 
this  is  an  expansion  ora  contraction  of  the  oval  with  the  change  in  the  Q  index.  It 
would  be  necessary  (but  rather  difficult)  to  separate  these  changes  for  determining 
the  exact  relations  in  the  corresponding  values  of  probabilities  in  Tables  5  and  6. 

Two  factors  still  missing  in  the  above  discussion  are:  (1)  a  varying  time  widtii 
of  the  ERS  test  area  at  different  CG  latitudes;  and  (2)  the  time  (CG1.T)  dependence 
of  probabilities  due  to  systematic  diurnal  motion  of  the  auroral  oval. 

For  determining  the  probability  of  auroral  disturbances  over  a  wide  longitude/ 
time  zone  covered  by  a  test  area  (like  one  in  Figure  7),  one  has  to  consider  the 
length  (w'idth)  of  the  window  subtended  by  the  sector  at  each  corrected  geomagnetic 
latitude.  Effectively  it  is  the  angle  subtended  by  each  arc  of  a  CG  Latitude  at  the 
CG  pole.  The  measured  angles  at  various  CG  latitudes  in  the  sector  of  Figure  7, 
are  presented  in  Figure  8.  In  Figure  8,  the  abscissa  is  the  CG  latitude  and  the 
ordinate  on  the  left-hand  side  is  the  angle  subtended  by  the  arc  of  a  CG  latitude  at 
the  CG  pole.  The  ordinate  on  the  right-hand  side  of  the  figure  presents  the  corres¬ 
ponding  durations  in  hours.  From  Figure  8  it  is  seen  that  the  time  width  of  the 
zone  is  less  than  15  min  for  CG  latitudes  below  4ii°  and  for  CG  latitudes  above  78°. 
The  time  width  of  the  sector  reaches  a  maximum  of  250  min  for  latitudes  59°  to  61°. 


Figure  8.  Time  Coverage  at 
Different  Corrected  Geomagnetic 
Latitudes  of  the  ERS  Test  Area 
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For  the  diurnal  motion  of  the  auroral  oval,  the  latitude  dependence  on  the  Q 
index  and  on  the  CGLT  is  computed  from  Starkov  Eq.  (1)  and  is  shown  in  Figure  9. 

In  Figure  9  the  abscissa  is  the  corrected  geomagnetic  local  time  (CGLT)  and  the 
ordinate  presents  the  equatorward  boundary  of  the  oval.  The  curves  are  for  various 
values  of  the  Q  index.  From  the  figure  it  is  seen  that  ovals  with  Q=0,  1  do  not 
reach  the  latitude  A=65°.  Even  for  the  highest  Q=9,  the  oval  is  poleward  of  A=65° 
for  the  duration  0815  to  1815  CGLT.  Thus  even  in  the  case  of  the  strongest  mag¬ 
netic  activity  Q=9  with  a  continuous  duration  of  24  hr  the  oval  would  reach  A=65° 
only  for  14  hr  from  1815  to  0815  CGLT.  The  latitude  would  be  free  of  auroral 
disturbance  for  the  remaining  10  hours.  This  systematic  diurnal  effect  has  to  be 
incorporated,  so  that  the  probabilities  in  Table  6  would  provide  a  basis  for  opera¬ 
tional  use. 


CORRECTED  GEOMAGNETIC  LOCAL  TIME 

Figure  9.  The  Diurnal  Latitude  Dependence  of  the 
Auroral  Oval  for  Various  Q  Values 


For  operational  use.  one  would  seek  the  probability  distributions  at  various  CG 
latitudes  under  varying  assumptions  such  as:  (a)  What  is  the  probability  for  the 
auroral  oval  to  affect  a  specific  latitude  and  longitude?:  (b)  What  is  the  probability 
that  the  whole  sector  in  the  test  area  along  a  given  CG  latitude  is  affected  bv  the 
auroral  oval?;  or  (c)  What  is  the  probability  that  the  auroral  oval  would  affect  any 
section  in  the  test  area  along  a  given  CG  latitude? 

Toward  achieving  this  goal,  first  the  general  case  a'  of  an  auroral  oval 
affecting  a  given  latitude  and  a  given  longitude  will  be  considered  in  detail.  This 
will  be  followed  by  examples  of  two  selected  latitudes  for  answering  the  questions 


J 
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b'  and  c'.  One  would  then  be  able  to  follow/use  the  procedure  for  deriving  the 
probabilities  for  any  latitude  covered  by  the  ERS  test  area  for  satisfying  the  speci¬ 
fic  needs  like  those  for  'b  or  'o'. 

For  the  case  'a'  the  probabilities  from  corrected  Q  in  Table  2f  are  combined 
with  the  systematic  diurnal  motion  of  the  auroral  oval  seen  in  Figure  9.  These 
probabilities  are  presented  in  a  series  of  graphs  in  Figures  10-1  to  10-18  to  cover 
the  CG  latitude  range  of  58°  to  75°. 

Figure  10-8  is  for  latitude  A=65°.  A  reference  to  Figure  9  shows  that  a  mag¬ 
netic  disturbance  with  Q2  2  is  needed  for  the  auroral  oval  to  icach  this  latitude.  A 
reference  to  Table  2f  shows  that  the  highest  cumulative  probability  (column  2)  for 
Q  -  2  is  68  percent.  Thus  the  highest  probability  in  Figure  10-8,  the  height  of  the 
envelope  is  68  percent.  Starting  from  the  corrected  magnetic  local  noon,  the  oval 
first  reaches  this  latitude  at  1645  CGLT  and  stays  (or  goes  further  equatorward)  up 
to  0915  CGL.T.  During  the  period  0915  to  1645  CGLT  the  oval  is  poleward  of,  and 
does  not  reach  the  latitude  A=65°.  Also  it  is  recalled  that  in  Figure  9  for  the  highest 
value  Q=9  the  duration  of  the  oval  at  any  single  longitude  and  with  A=651’,  is  14  hr 
(1815  to  0815  CGLT).  In  Figure  10-8  the  maximum  duration  is  (1645  to  0915) 

16-1/2  hours.  This  width  is  based,  in  addition,  on  the  distribution  of  the  corrected 
Q  index  of  1974.  The  first  curve  on  the  left-hand  side  in  the  envelope  is  for  a 
minimum  continuous  duration  of  15  hr  for  the  disturbance  at  A=65°.  The  last  curve 
on  the  right-hand  side  inside  the  envelope  is  for  the  duration  of  disturbance  to  last 
at  least  for  1  hr  at  A=65°,  For  example  once  the  CGI.T  is  07  15  at  A=65°,  the 
longest  duration  for  this  latitude  to  be  affected,  cannot  be  more  than  2  hr  (from 
0715  to  0915  CGLT). 

Figure  10-15  presents  the  expected  probabilities  at  A=72°.  The  Starkov  Eq.  (1) 
shows  that  normal  .  the  oval  would  reach  this  latitude  (from  the  pole)  for  12  hr  if 
Q=0.  In  the  figure  the  100  percent  probability  is  from  1830  to  0800  CGLT,  for  a 
period  of  12-1/2  hours.  .  The  increased  duration  arises  from  the  fact  that  the  higher 
Q  values  push  the  oval  to  lower  latitudes  (<  77°)  in  the  night  sector,  in  the  figure, 
the  probability  never  drops  below  42  percent  at  this  latitude.  This  minimum  prob¬ 
ability  is  only  for  the  period  from  1100  to  1500  CGLT. 

In  the  above  discussion  the  results  are  based  on  the  oval  boundary  determined 
from  discrete  auroral  arcs  from  the  Starkov  equation.  An  invading  aircraft  could 
fly  unnoticed  under  the  auroral  oval,  manifesting  a  threat  over  the  ERS-area. 

From  the  use  of  graphs  in  Figure  10,  one  can  predetermine  the  levels  of  relative 
threat  at  various  CG  latitudes  with  time  (CGLT)  along  any  given  longitude  of  the 
ERS  test  area. 

For  answering  the  questions  ’b'  and  'c'  one  has  to  take  into  account  the  time 
widths  at  respective  CG  latitudes  of  the  ERS  test  area.  These  time  widths  have 
already  been  presented  in  Figure  8. 
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When  one  considers  the  threat  over  the  entire  length  along  a  given  CG  latitude 
of  the  test  area,  the  time  intervals  are  limited  such  that  the  auroral  oval  simul¬ 
taneously  affects  both  the  eastern  and  the  western  edges  of  the  test  area  along  that 
latitude.  This  additional  constraint  has  to  be  introduced  in  the  graphs  of  Figure  10. 

Let  us  consider  the  situation  at  A  =  65°.  For  the  sake  of  convenience  Figure 
10-8  is  reproduced  as  Figure  llal.  It  is  recalled  that  the  graphs  in  Figure  10 
referred  to  only  single  longitude.  Referring  to  the  ERS  test  area  in  Figure  6,  the 
results  can  be  applied  to  the  corrected  geomagnetic  local  time  of  the  western 
longitude  of  the  ERS  test  areas,  as  well  as  to  its  eastern  longitude.  Though  the 
CG  Local  times  are  the  same  for  both  the  western  and  the  eastern  edges  of  the  area, 
the  respective  local  times  are  separated  from  each  other  by  250  min  (from  Figure  8). 

When  the  condition  of  simultaneous  effect  at  both  the  ends  (east-west)  is  imposed 
the  time  width  of  the  probability  curve  would  reduce  (symmetrically)  by  that  amount. 

The  corresponding  curves  are  presented  in  Figure  1  lb  1 .  In  Figure  1 1  b  1  the  height 
of  the  probability  curve,  being  controlled  by  the  distribution  of  the  corrected  Q 
index,  is  unchanged.  The  width  of  the  peak  is  reduced  to  45  min  (0045  to  0130) 
compared  to  that  of  285  min  (2245  to  0300)  in  Figure  llal.  The  width  at  the  base 
90  percent  probability)  is  reduced  correspondingly  from  16-1/2  in  Figure  llal  to 
12-1/2  hr  in  Figure  1  lb  1. 

Figures  lla2  and  llb2  present  results  for  latitude  A- 72°.  The  results  are 
similar  to  those  for  A=65 0  in  Figures  llal  and  1 1  b  1 .  In  both  Figures  lla2and 
1  lb2,  the  minimum  probability  for  auroral  disturbance  is  42  percent.  The  mere 
motion  of  the  Starkov  oval  indicates  that  the  duration  should  be  12  hr  for  100  per¬ 
cent  probability  due  to  auroral  disturbance  at  A=72°.  It  is  already  shown  that  the 
Q  index  distribution  increased  this  duration  (Figure  1  1  b  1 )  to  13-1/2  hours.  Refer¬ 
ring  to  Figure  8,  the  time  width  of  the  ERS  test  area  at  A=72°  is  iBO  minutes.  This 
reduces  the  width  o'f  simultaneous  effect  over  the  test  area  at  A=7  2°  by  3  hr  from 
Figure  lln2to  figure  llb2. 

In  answering  question  'o'  it  follows  that  one  would  have  to  extend  the  time  width 
such  that,  even  if  one  edge  (eastern/western)  is  clear  of  any  auroral  disturbances, 
the  possibility  that  the  remaining  edge  of  the  test  area  could  be  affected,  has  to 
be  considered.  Thus  the  widths  of  the  graphs  in  Figure  10  have  to  lie  increased  to 
take  into  account  the  time  width  of  the  ERS  test  area  at  various  latitudes  already 
shown  in  Figure  fl.  The  corresponding  results  are  presented  in  Figures  llcl  and 
1  lc2  for  latitudes  !> fa °  and  72°  respectively.  In  contrast  to  the  reduction  in  widths 
from  Figures  llal  md  11, u2  to  Figures  llbl  and  llb2,  one  now  sees  ar,  increase 
in  the  widths  bv  the  same  amounts  in  tin-  respective  Figures  llcl  and  1L-2  (from 
f  lgures  llal  and  1  I a2). 

rims  one  can  select  and/or  derive  a  set  of  graphs  from  f  igure  10,  to  satisfy 
the  specific  needs  different  from  those  of  the  operational  use  of  the  graphs  in 
figure  10. 
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Figure  llbl.  Probability  That  the  Entire  ERS  Test  Figure  llb2.  Similar  to  That  of  F'igure  llbl  for  A=72 

Area  at  A=65°  Would  be  Affected  by  the  Auroral  Oval 


4.  CONCLUSIONS 


The  presence  of  one  or  more  phenomena  such  as  discrete  aurora,  continuous 
(diffuse)  aurora,  sporadic  E  layer  E  ,  auroral  E  layer  E  ,  F„  layer,  F  layer 
irregularity  zone  (FLIZ),  the  ionospheric  trough,  could  affect  the  test  area.  De¬ 
pending  on  the  factors  such  as  the  locations  of  the  transmitter  and  the  receiver  of 
the  radar  system,  the  emission  power  of  the  transmitter,  the  threshold  signal 
detectability  of  the  receiver,  one  has  to  determine  as  to  which  of  the  above 
phenomena,  the  radar  system  is  most  susceptible. 

The  statistical  study  of  the  Q  index  data  in  1974  from  Sodankyla,  in  conjunction 
with  the  ERS  test  area,  shows  that  only  about  one-fourth  of  the  poleward  sector  of 
the  area  is  affected  by  auroral  disturbances  for  more  than  50  percent  of  the  time. 
The  middle  one  -third  portion  is  affected  for  less  than  25  percent  of  the  time  and 
the  bottom  third  portion  is  free  of  any  auroral  disturbances  at  all  times. 

In  all  the  above  discussion,  the  location  of  the  equatorward  boundary  of  dis¬ 
crete  auroral  arcs  from  the  Starkov  equation,  was  the  criterion  affecting  the 
regions  in  'he  test  area.  For  example,  if  instead,  the  radar  system  is  susceptible 
to  diffuse  aurora,  in  preference  to  discrete  auroral  arcs,  another  study  has  shown, 
that  the  diffuse  aurora  statistically  lie  3°  equatorwards  of  the  discrete  arcs.  Thus 
the  probabilities  presented  in  Tables  5  and  6  at  various  latitudes  have  to  be  moved 
3°  equatorwards.  Thus  the  present  probability  at  75°  CGL  from  discrete  auroral 
arcs  has  to  be  used  for  7  2°  CGL  in  case  of  diffuse  aurora. 

This  report  has  developed  a  technique  to  determine  the  effect  of  the  variation 
of  the  Q  index  measured  from  a  single  auroral  station  on  the  radio  wave  surveil¬ 
lance  in  the  coverage  area  of  a  OTH  backscatter  radar  system.  To  facilitate  this 
demonstration  a  number  of  illustrations  and  tables  have  been  presented.  This  will 
ease  the  application  of  this  technique  to  any  other  test  area  covered  by  any  radar 
system. 

The  first  section  shows  that  the  Q  index  measured  from  a  single  auroral  station 
exhibits  a  diurnal  variation  due  to  the  movement  of  the  earth  under  an  auroral 
oval.  This  systematic  diurnal  variation  can  be  eliminated  effectively  by  deter¬ 
mining  a  set  of  empirical  factors.  The  corrected  data  set  from  a  single  station 
may  then  be  used  reliably  for  computing  the  durations  as  well  as  the  levels  of 
magnetic  (and  auroral)  disturbances. 

The  second  section  combines  the  location  of  the  auroral  oval  at  various  (.'(! 
latitudes  expected  from  the  Starkov  equation,  with  the  corrected  durations  and 
corrected  levels  of  the  magnetic  activity  index  <.).  for  obtaining  file  probabilities 
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for  the  iii  st  i.  rli  am  es  global  base.  These  probabilities  are  ..mailer  for  a  loca¬ 

tion  it  its.  single  selected  longitude.  1  he  lahu  s  for  probabilities  for  a  test  area 
’.\ith  .1  ilefipitr  b  Tiglt  aim  il  i  It  ;  lie  between  these  two  eXtiemes  (of  cables  5  and  6). 

In  the  )TH-l!  lilt-*  coverage  •  r<  1  the  window  varies  with  the  '01  atiuns  along 
the  CC.  l  ititufles  t  s  sh>*v.ti  in  1  ig.i  j  e  l  os  one  rational  .s.ige  the  combined  effect 
of  the  time  width  of  tr.e  K Hr  U  st  .  ,t  •.  .riu..n  latitude..,  the  svstemutie  diurnal 
motion  ef  tne  oval  .1:1  1  ri<  t 1  11  st :  ib..* cm  .  the  eo r ta-cted  index  ■  1 : * 7  4  •  is  *  1  ken 

into  account.  The  final  results  for  pricti*  ’  is*  are  presented  is  ■.  s*  o' 

graphs  (figures  10-1  through  I  0  -  1  rf  >  for  variou.  latitudes. 

This  detailed  step-bv-step  presen'  t»ion  si. on!  ;  enable  one  to  use  the  v,*  index 
data  over  an”  period  lor  determining  the  expo.  *•  u  probabilities  of  am  1  .!  dis  • 
turbances  for  different  h  els  md  dur  itions  of  ’.umnetii  utility  ''or  at-  ted  .ire.* 
to  be  under  surveillance. 
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